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Abstract: We point out that there exist two different formulations of the Heavy 
Quark Effective Theory (HQET). The one formulation of HQET was mostly de- 
veloped at Harvard and involves the use of the equation of motion to eliminate 
the small components of the heavy quark field. The second formulation, developed 
in Mainz, involves a series of Foldy-Wouthuysen-type field transformations which 
diagonalizes the heavy quark Lagrangian in terms of an effective quark and anti- 
quark sector. Starting at 0(l/mQ) the two formulations are different in that their 
effective Lagrangians, their effective currents, and their effective wave functions dif- 
fer. However, when these three differences are properly taken into account, the two 
alternative formulations lead to identical transition or S-matrix elements. This is 
demonstrated in an explicit example at ©(l/mg). We point to an essential diffi- 
culty of the Harvard HQET in that the Harvard effective fields are not properly 
normalized starting at order 0{l/rn^). We provide explicit higher order expressions 
for the effective fields and the Lagrangian in the Mainz approach, and write down 
an ©(l/mg) nonabelian version of the Pauli equation for the heavy quark effective 
field. 
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Let us get straight into the heart of the matter of what we want to discuss in 
this talk by stating that there exist two different formulations of the Heavy Quark 
Effective Theory (HQET) which will be referred to as the Harvard HQET Q and 
the Mainz HQET 0, |. 

The HQET that is most widely in use was mostly developed at Harvard 
It is for this reason that we shall refer to this version as the Harvard HQET. In 
the Harvard approach one first extracts the mass phase from the heavy quark field 
iPq{x) and then splits up the residual field into its "large component" and "small 
component" pieces h{x) and H{x), respectively. Accordingly one writes 

^PQix) = e-'"'^^■''ihix)+H{x)) (1) 

where 

j^h{x) = h{x) (2) 

and 

]^H{x) = -H{x) (3) 

and where is the four- velocity of the heavy quark, = p^/mq. Unfortunately 
the nomenclature in terms of the "large" component field h{x) and the "small" 
component field H{x) has been somewhat tangled up in the course of developing 
the Harvard HQET. In the following we shall drop reference to the x-dependence in 
the fields. 

In the next step one takes "one-half of the full equation of motion of QCD, i.e. 

{i]p-mQ)^Q = Q (4) 

by applying the small component projector (1 — ;^)/2 to Eq. (^): 

[iv D + 2mQ)H = i^^h (5) 

In Eq. (|^) we have introduced the four-transverse component of the covariant deriva- 
tive 

= D^ -V Dvf" (6) 

where the four-transversality is defined w.r.t the velocity v^. Eq. (^) can then be 
inverted to obtain 

H = 7 ^ipih (7) 

The small component field H can be seen to be related to the large component field 
h through a geometric series in the inverse power of the heavy mass. The Harvard 
approach consists in eliminating the small component field out of the theory via Eq. 
(0). Before we proceed any further there are two asides that we want to embroider 
on. 

First, there is the concept of the four-transversality used in Eqs. (|^) and (0). 
Note that any four-vector can be split into its components transverse and parallel 
to a given four-velocity according to 

= a^A^ + af := {a^" - v av^) + {v ■ av^") (8) 
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The concept of four-transversality is an important concept in the formahsm of HQET 
as one is frequently referring to rest-frame objects where = (1, 0, 0, 0) and where 
four-transverse vectors reduce to pure three-vectors a±'^ a. This concerns the 
covariant derivative in Eq. (|^) as well as the relative momenta ki'^ = k'^ — k-vv'^ and 
spin operators 'ji^ = jt^ — ■^yi^ that are needed in the construction of the spin wave- 
functions of HQET 1^. Explicit reference to the /onr-transversality is important in 
order not to get confused with the usual notion of transversality which refers to a 
three-transversality , i.e. a = a — a ■ kk {k is & three-momentum of unit length). 
In fact, the concept of four-transversality should be quite familiar from QED when 
one is considering virtual photon exchange. In the rest frame of the virtual photon 
q = the conserved vector current amplitude = (6| J^|a), with g'^ = 0, reduces 
to a three-vector corresponding to the spin-one nature of the virtual photon. When 
expanding the conserved amplitude in terms of covariants in a general frame this 
has to be done in terms of four-transverse covariants e.g. p_L^ = — ^g^, where 
the four-transversality is defined w.r.t. the momentum of the photon. 

The second aside concerns an alternative derivation of the Harvard HQET that 
was proposed in 0. The authors of employ functional integration techniques to 
integrate out the small component field H from the functional action. They then 
arrive at the same relation Eq. (|^). That the equation of motion and functional inte- 
gration derivations of HQET are entirely equivalent may be appreciated by consider- 
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ing the following simple example using the classical Lagrangian L{x, y) = ~ \ ~\~ U^- 
If one wishes one may view x as being related to the small component field H and 
y as been related to the large component field h. Next consider the classical action 
integral 

/ dy r dxe^^'^^y^ = fdy r rfxe-^"'+^" (9) 

J J—oo J J —oo 

The X- integration in (^) can be done by the usual "completion of the square" trick, 
i.e. by writing — + yx = —^^^^ — |- ^ One then arrives at 

/ dy r dxe^'-'^'y^ = V2^ f dye"^ (10) 

J J — oo J 

After having integrated out the x-degree of freedom the new Lagrangian reads 
L{y) = y^/2. The crucial observation is that the same Lagrangian is obtained by 
using the equation of motion obtained by variation of the L[x, y) w.r.t. the x-degree 
of freedom, i.e. dL{x, y)/dx = 0. This then gives the equation of motion —x + y = 0. 
Substituting for x in the original Lagrangian L{x,y) one obtains the same new La- 
grangian L{y) as before. We have chosen this simple illustration in order to make 
the point that the equation of motion and functional integration derivation of the 
Harvard HQET are entirely equivalent. This can be appreciated without having to 
go into the technicalities of the full functional integration approach presented in . 
Returning to Eq. (0) one writes 

H=^ ^ \ ^ m^h = ^(l-'^+ ^'''-'l^\ ..\rV)^h (11) 
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As mentioned before one then ehminates the small component field H by substituting 
Eq. (0) into the effective Lagrangian and current expressions. This defines the 
Harvard HQET in terms of a l/mq geometric type series expansion of the effective 
Lagrangian and current in terms of the large component field h only. 

In the Mainz HQET [§, ^ one uses a series of exponential Foldy-Wouthuysen- 
type field transformations which yield an exponential type series expansion for the 
effective Lagrangian and the effective current. 

At this point we want to be rather suggestive and write down geometric and 
exponential series expansions for a given small parameter a. One has : 

geometric series: = 1 + m — ■ • • (12) 

1 — ia 

exponential series: e*"^ = 1 + ia — -a^ ■ ■ ■ (13) 

The series have been arranged such that they coincide in the first order term. Al- 
though the two series would be an oversimplified representation of the two effective 
theories it still makes the point correctly: The Harvard and Mainz HQET's are 
different in that their effective Lagrangians and currents are different starting at 
0{l/mQ). However, despite the fact that the respective effective Lagrangians and 
currents are different, one expects the physics i.e. physical on shell matrix elements 
to be the same in both theories, if they are calculated correctly. That these issues 
are of immediate practical concern is being evidenced by the fact that l/rriQ and 
even higher order corrections are presently being disscused in the literature 0. 

As a next step we want to explain in somewhat more detail how one derives the 
Mainz HQET. To start with let us briefly retrace the physics steps that lead one 
from the full QCD Lagrangian to the static heavy quark Lagrangian. This will first 
be done in a completely heuristic manner. Consider again the full QCD Lagrangian 

C = ^Q{zP-mQ)^Q (14) 

The static approximation consists in neglecting the three-derivative in Eq. ([T^ 
relative to rriQ according to the expansion 



E=^ml + P = mQ (^l + £^ + ...^ (15) 

This can be achieved by an appropriate field redefinition ipq ip'q. This will then 
lead to 

Step I: = ijj'Qii-foDo - mQ)^)^ (16) 

In the next step one shifts the energy scale E —* E' = E — niq which can again be 
achieved by an appropriate field redefinition ip'q hq. The Lagrangian now reads 

Step II: = hq{i-foDo)hq (17) 

If one chooses to work in terms of the heavy quark field h^~^^ only (with h^^^ = 
|(1 + 7o)^q) one then recovers the static Eichten-Hill Lagrangian C = TS^^HD^h^'^'' 
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The step-wise reduction of the QCD Lagrangian (JT^) to the final form (|T^) can be 
achieved by a series of Foldy-Wouthuysen-type field redefinitions which eventually 
yields the leading term result ( ]T7D as well as all higher dimension operators in the 
l/mg expansion 0. The first transformation is (j = 1,2,3) 

Step I: ^ e*^^^^/2'"«^^ 

V'Q ^'Qe-^^^^j/^"« (18) 

where the arrow on the derivative indicates in which direction the derivative acts. 
The heavy quark Lagrangian now becomes 

oo / \ \^ 



giving a form of the Lagrangian which makes explicit the mass perturbations. 

The second transformation that removes the heavy mass dependence in the first 
term of (|19D can be seen to be given by 

Step II: tfjQ ^ e-'^^Q^'^hg 

/iQe™«^»* (20) 

However, in order not to bring in further numerator mass terms through the 
transformation (^) one needs to first block-diagonalize the higher order operators 
Ok appearing in ([T9|) . The block-diagonalization has to be done w.r.t. the upper 
and lower component of the heavy quark fields ^-y^^Q and ^—^ipq, respectively. 
This can be achieved order by order by splitting the operator Ok into two pieces O^ 
and O'^ that commute and anticommute with 70, respectively. That is, one writes 

Ok = 01 + 01 (21) 

where 

Or = \{Ok±ioOklo) (22) 

and 7o] = and {C^, 70} = 0. The anticommuting operator O^ is then removed 
from that order of {l/rnqY by a further exponential type transformation. It is 
literally shifted down to become part of the higher dimension operator Ok+i- 

Note that one never introduces any further implicit time-derivatives do through 
the field redefinitions in addition to the time- derivative in the lowest order term of 
the final Lagrangian (see (|T^). Technically speaking, explicit higher order time- 
derivative terms do appear through the above field redefinitions. However, these 
higher order time-derivatives always come in as commutators that are related to the 
field strengh tensor (see Eq.(31))0, This will become important when we discuss 
the Pauli equation later on. 

Up to now we have remained in the rest frame = (1, 0, 0, 0) of the heavy quark 
(or antiquark) in order to stay as close as possible to the heuristic considerations that 
led us to the static Lagrangian (|T^. The field redefinitions that lead to the HQET 



4 



Lagrangian can in fact be done in any moving frame with 7^ by effecting the 
replacements 



70^0 ^ :=^|| 
--fjDj ^ ]/) - ■ D := ]/)^ (23) 

In the moving frame the heavy quark effective Lagrangian is then given by 
^hIet = hQ^v ■ D + ^i-D' + {v . Df 

+ '-gcj^,F^^ - tg^.^F^'v,) + ■ ■ ■}hQ (24) 

where F^'^ is the field strength tensor F^" = [D^, D^] {D^ = — igA^^ and A = 
A^-t"-) and where we have kept terms up to 0{l/mQ) only. We chose to label the 
Mainz effective theory by the initials of the authors of 0. Note that the HQET 
Lagrangian contains both heavy quark fields h^^^ and heavy antiquark fields h^~^ 
where ^/i'-^^ = i/i^^-' and Hq = h^^^ + h^~\ Remember that this is different in the 
Harvard approach where the construction is done either in the quark sector or in 
the antiquark sector by adjusting the sign of the mass phase in Eq.(l). The quark 
and antiquark sector in the Mainz HQET are completely decoupled (at any order!) 
which is guaranteed by the very procedure of deriving the HQET Lagrangian. To be 
sure one can easily convince oneself that the nondiagonal contributions h^^'^ . . . /i*^"^ 
and /z'-"^ . . . h^'^^ induced by the first order terms a^^F'^'^ and 'j^j^F^'^Vi, in Eq.(p4D 
cancel. On the other hand, the j^i^F^'^Vi, contribution vanishes for the diagonal 
contributions h^~^^ ■ ■ ■ h^^^ and h^~^ ■ ■ ■ h^~^ and one thereby recovers the 0{l/mQ) 
Harvard HQET Lagrangian as promised before (see Eq. (|26|) ). 

As noted before the method of deriving the Mainz HQET proceeds step-wise 
order by order. This is an iterative procedure which lends itself to computer imple- 
mentation. In fact we have written an efficient program in Mainz that achieves just 
this. For the fun of it the program has been run up to 0{l/m^) |0. As a sample 
result we list the Mainz HQET Lagrangian for the heavy quark field h^'^^ up to order 
O^I/ttiq). One has 



riCT _ T(+) 
1 



2 771 I 771 771 2 J 



^^ll^^^ll^^ +^||^^^^|| +^^^11^^^ 

^ ^ 'im m m \ ^ m "im "in i'^77i77i77i2i 



16m} 
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32mQ 



10„ o„ „ „ 4 



+ + 2^^^)] (25) 

As mentioned before, the time-derivative terms appearing at ©(l/mg) and higher 
can all be rewritten in terms of the field strength tensor 0. 

The ensuing discussion will be in terms of the heavy quark field h^'^^ only and, for 
the sake of convenience, we shall omit the label (+) on the heavy quark field in the 
following. In order to pinpoint the differences in the Harvard and Mainz HQET's let 
us consider the heavy quark effective Lagrangians and the flavour-conserving heavy 
quark effective currents up to 0{l/m'Q). As mentioned before it is at this order 
where they start to differ from one another. For the Harvard HQET one has 

qTt' = h{^v ■ D - + ^P^v ■ DPJh (26) 



iT^vD^^+^^T^^^+^^vDT^)}h (27) 



Am, 

The 0{l/mQ) difference terms for the Mainz HQET are given by |10 



^HQET = ^ZTet' + [-^P^'v ■D--vDP^')h (28) 



j;^,HQET = J" HQET ' (^I^/e, + ,P^'j h (29) 

There are a number of observations we want to make about the difference terms in 
Eqs. (H) and (|2|). 

Let us first rewrite the difference term in the effective Lagrangian (PH). One has 



- -D-^v DP^' = -P^v ■ DP^ - ^PAA, v-D] + ^[P^,v D]P^ (30) 
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When one looks at the time derivative term v ■ D in (pOD (or alternatively Dq in 
the rest frame) one can see that the first term on the r.h.s. of ( ^O]) cancels the time 
derivative term in the Harvard Lagrangian The remaining two terms on the 

right hand side of ( ]28|) do not contain true time-derivatives since the commutator 
can be expressed in terms of the field strength tensor via 



(31) 



The difference in the Mainz and Harvard Lagrangians lies in the fact that the time- 
derivative terms only appear at leading order in the Mainz Lagrangian but to all 
orders in the Harvard Lagrangian. This is exemplified at second order in the above 
example. As a consequence of this one can therefore quite easily write down Pauli 
equations to any desired order in the Mainz approach by using the equation of 
motion for the effective heavy quark field h in the rest frame of the heavy quark. 
For example, this has been done to (9(l/mQ) in The result reads 



ot \ 2m 2m 



■a ■ TotE^t" 



im^ 



im^ 



im-' 



[{P - gA'^ey - ga ■ BH- 

EH" ■ EH^ + -fabca ■ {E" X E^)e 
2 



^[dwE- + Uc^ ■ E')t'^ 
a ■ X (P - gA^t^) 



am 
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Am^ 



h 



(32) 



where E and B are the electric and magnetic colour fields. To our knowledge the 
Pauli equation has not been derived to this order before, let alone in the non-Abelian 
case. We hope that we have convinced the reader by now that the Mainz approach 
to HQET allows one to do so very efficiently. 

As the next topic we want to discuss the calculation of physical matrix elements 
using the two different formulations of HQET. We again concentrate on the 0{1/ mq) 
contributions where the two formulations start to be different. 

and (|29|) in a rather sug- 



We have arranged the difference terms in eq. (^ 



gestive manner by placing them one below another: except for the v ■ D factors 
in the Lagrangian the difference terms in the Lagrangian and in the current look 
quite similar to one another. In fact one finds that the differences in the effective 
Lagrangians "almost" cancel in physical transition matrix elements. They cancel 
except for some 0{l/mQ) contributions which can be absorbed in the definition of 
the 0{l/mQ) HQET wave functions or the interpolating field associated with the 
HQET state. 

Let us explicitly demonstrate this cancelation in the fiavour-conserving case by 
using the Feynman diagram language. One of the second order contributions of the 
Lagrangian difference term is drawn in Fig. la. for b ^ b transitions. The v ■ D term 
from the "nonlocal" Lagrangian insertion can be seen to cancel against the inverse 
propagator {v ■ D)^^ adjoining it. In the coordinate space language the two vertices 
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in Fig. la. are contracted to one point due to the relation {iv ■ D) S{x, y) = S{x — y). 
One thus remains with the effective "local" insertion drawn in Fig. lb. This will be 
canceled by the truly "local" l/rn^ insertions from the effective current difference 
(p9|) which is not drawn here. Nonlocal insertions in which v ■ D operates on the 
outer heavy quark wave function vanish due to the equation of motion. There is one 
exception to this which is represented by the Feynman diagram drawn in Fig.lc. The 
v-D to the left of the nonlocal insertion is now canceled by the inverse propagator to 
the right of the local insertion and one remains with a 0{l/m^) difference even for 
physical transitions. This difference corresponds to a 0{l/mQ) renormalization of 
the HQET wave function or, more precisely, the interpolating field that is associated 
with the HQET state. The upshot of the analysis is that the Mainz and Harvard 
formulations of HQET lead to identical physical matrix elements as long as one 
takes into account the differences in the 0{l/mQ) wave functions as drawn in Fig. 2 
in addition to the differences in the effective Lagrangians and effective currents. The 



same conclusion was reached in |]T0| using functional differentiation techniques. 

However, in spite of the fact that the predictions for physical matrix elements 
are identical, there exists an essential difference between the two HQETs, in that 
the Harvard HQET lacks a correct normalization of the heavy quark field. We will 
illustrate the problem in the context of QCD with unconfined heavy quarks. When 
evaluating S-matrix elements for scattering processes by using the usual machinery 
of the reduction formalism, an assumption is implicitly made about the value of 
matrix elements like 

{Q\hixmp,s)) = ^^V^«(^Q^'^)e-''■^ (33) 

giving the wave-function of the ingoing or outgoing heavy quarks. This prejudice 
has its roots in the fact that such a relation holds indeed true (by construction), if 
the heavy quark field in the effective theory h{x) is replaced by the QCD field Q{x). 
However, this relation is not true if h{x) is the heavy quark field in the Harvard 
HQET. The right-hand side has to be multiplied with a correction factor of the 
form 

' + 8^ + --- ^^^^ 

which deviates from unity by an amount proportional to the residual momentum 
p of the one-quark state \Q{p,s)). It is clear that the omission of this correction 
factor, such as it is done in the naive application of the Harvard HQET, will possibly 
lead to a wrong answer. 

On the other hand, relation (33) is valid to any order in the Mainz HQET. This 
has been shown in by proving that the heavy quark field in the Mainz effective 
theory is connected to the QCD field by a unitary transformation. The fact that 
the Harvard HQET field has a different normalization can be appreciated by noting 
that the two are related by 

/I^^(X) = f 1 + h^-rvard ^35) 

V 8?77,r) / 
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Another way of looking at the issue of the normahzation of the heavy quark 
field is to examine the form of the heavy quark propagator in the vicinity of the 
one-particle pole. The requirement of correct normalization can be expressed by 
saying that the residue of this pole should be unity. That these two points of view 
are related can be seen from the form of the Kallen-Lehmann representation for the 



heavy quark propagator ^1 



Spip) = ^ / dM I ■ A 36 

2 Jo i{v ■ p — + le) 

where p{M'^) is defined by a sum over all possible intermediate states 

p(/)5„^ = {2nfY.S^'\Pn-p){0\ham'n^){n\hm\0) . (37) 

n 

If one then makes use of Eq.(33), the contribution of the one-quark intermediate 
state to the spectral function p(M^) can be seen to be as follows 

p{M^) = 6{M^), (38) 

which gives a free field-like behaviour of the heavy quark propagator in the neigh- 
bourhood of the one-particle pole. On the other hand, in the Harvard HQET the 
residue of the propagator at the one-particle pole takes a different value (which is 
even a function of the residual momentum of the heavy quark) . This has been shown 
through explicit calculation in a recent paper by A. Das [|T^ . 

The lack of maintaining the correct normalization of the field or wave functions 
by eliminating the "small" component via the "equation of motion" approach has 
long been known in the Abelian context of QED, where these issues were studied in 



connection with the Foldy-Wouthuysen transformation |T^, |13| [IJ, |15|. Eliminating 
the "small" component naively leads to a nonhermitean Hamiltonian. In the case of 
an electron interacting with an external electric field, the lowest order manifestation 
of the non-hermiticity is an imaginary electric dipole moment. All these issues were 
very nicely discussed in the paper by A. Das [T^. 
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Figure Captions 

Fig.l 0{l/m^) contributions from the effective Lagrangian to the flavour-conserving 
H},{v) — > Hb{v') current transition, a) and b) show how the nonlocal contribution a) 
becomes a local current contribution b) through the contraction of two interaction 
points. In c) we show an external line contribution which can be absorbed into the 
definition of the HQET wave function. 
Fig. 2 0(l/mn) difference of Harvard and Mainz HQET wave functions. 
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